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ABSTRACT 

Using archival XMM-Newton observations of the diffuse and unresolved compo- 
nents emission in the inner disc of M33 we exclude the possible contribution from 
narrow line emission in the energy range 0.5 5 keV more intense than 10~ 6 -i- 10~ 5 
erg/s. Under the hypothesis that sterile neutrinos constitute the majority of the dark 
matter in M33, we use this result in order to put constraints on their parameter space 
in the 1 10 keV mass range. 
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1 INTRODUCTION 

One of the main problems of modern astrophysics and cos- 
mology is the unknown nature of the dark matter (DM). 
Despite the fact that it represents the 80% of the matter 
content of the Universe (Komatsu et al. |2011 1, it has so far 



escaped any unambiguous direct detection attempt in under- 
ground laboratories. This has given rise to the development 
of several detection techniques aimed at finding signatures 
of DM annihilation or decay in the spectrum of cosmic rays 
( Porter et al. |2011[ ). 

On the theoretical side, to solve this puzzle a wide 
range of possible solutions assuming new physics has been 
proposed. Standard particles are in fact unable to explain 
the evidence for the DM, while e.g., Supersymmetric and 
Kaluza-Klein extension of the Standard Model are known 
to provide a rich set of possible DM candidates in terms of 
weakly interacting massive particles. They have represented 
the standard assumption about the nature of DM for a long 
time, and even nowadays they constitute the preferred and 
most studied candidates in the literature. Nonetheless, many 
alternative solutions have been proposed. (See |Feng ||2010| 
for an up to date review.) 

Among these, a particularly interesting scenario as- 
sumes that DM is instead made up by 0(1) keV sterile 
neutrinos ( |Dodelson fc Widrow |1994| ). The well established 
discovery of neutrino masses points toward the possible exis- 
tence of gauge singlet fermions responsible for the neutrino 
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mass generation through a seesaw mechanism. The corre- 
sponding masses can range over all known mass scales, but 
if some of the singlet fermions are light, then sterile neutri- 
nos emerge in the low-energy effective limit of the theory. 
Such particles are predicted to decay into a lighter neutrino 
and a photon ( Pal fc Wolfenstein || 1982 I , producing a nar- 
row line in the X-ray spectrum at an energy equal to half the 
mass of the decaying neutrino. Such emission is expected to 



affect the formation of the first stars ( Biermann & Kusenko 



2006) and to be directly detectable by existing and future 
X-ray telescopes. Sterile neutrinos could also explain the ob- 
served velocities of pulsars thanks to their anisotropic emis- 
sion from a cooling neutron star born during the explosion 
of a supernova. (See Kusenko |2009 for a dedicated review.) 



Searches for sterile neutrinos have been conducted both 
in dwarf and giant galaxies, using all available X-ray obser- 
vatories. For instance |Loewenstein et al. | ( |200"9[ ) fail to find 
signatures of sterile neutrino decay, using the Suzaku ob- 
servations of the Ursa Minor dwarf spheroidal galaxy, and 
placed constraints on the existence of sterile neutrinos with 
given parameters. However, using the Chandra X-ray obser- 



vatory, Loewenstein & Kusenko (2010) reported the detec- 



tion of a narrow emission feature with energy of 2.5 keV 
in the dwarf spheroidal (dSph) Willman 1, consistent with 
emission line from sterile neutrino radiative decay. However 
using the same data |Nieto &: Mirabal | ( |2010| ), claim that 
all lines in the X-ray spectrum can be explained by resid- 
ual background emission and spurious features arising from 
an incorrect modeling of the latter. On the other hand[Bt>| 



yarsky et al. (2010) analyze a combined sample of XMM 
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Figure 1. The X-ray spectrum of the M33 inner region, from the 
PN camera; the best-fit model discussed in the text is overlayed, 
together with the individual components of the model. Note that 
only the two thermal components (those peaking at ~ 0.6 and 
~ 0.9 keV) are free to vary, while all the others (due to the bright 
source contribution) are fixed to the residual flux expected to leak 
outside the masked region. 



observations of M31, Willman 1 and Fornax dSph as well 
as Chandra data for Sculptor dSph, finding no evidence of 
a spectral feature at 2.5 keV, and confirming the previous 
findings of Watson et al. 2006| (See, however, Kusenko & 



Loewenstein 2010 for a discussion of systematics.) Mira- 



bal (20101 further used archival Swift XRT observation of 



the Milky Way satellite Segue 1, an ultra-faint dwarf galaxy 



with extreme mass-to-light ratio (M/L > 1300, see Simon 
et al. |201l] and |Geha et al. |2009| ). They find no evidence of 
emission lines due to neutrino decay either and put upper 
limits on its parameter space. 

In this work we make use of archival XMM~Newton ob- 
servation of the diffuse and unresolved components emission 
of the inner disc of M33 ( |Owen fc Warwick ||2010| >. We set 
constraints on the sterile neutrino parameter space and com- 
pare them with analogous results previously found in the 
literature. 

The paper is organized as follows: in section 2 we sum- 
marize the data set used in our analysis and the results we 
get; in section 3 we explore the consequences of our null de- 
tection in terms of bounds on the sterile neutrino parameter 
space; in section 4 we set our conclusions. 



2 X-RAY DATA ANALYSIS 

The X-ray data used in this work are those already dis- 



cussed extensively by Owen & Warwick ( 2010 1, and summa- 



rized here for completeness. The original observations were 
part of the |Pietsch et al. | ( |2004| > survey of M33. Among 
the entire dataset covering M33, we used only the obser- 
vations #0102640101, #0102642301 and #0141980801 from 
the EPIC-PN detector, i.e. the same used by |Owen fe War^| 
(20101 for their spectral analysis (see their Section 



wick 



2.3) and covering the central part of the M33 disk out to a 
distance of ~ 3.5 kpc, with observing times ranging from 7 
to 10 ks. The details of the X-ray data reduction (tempo- 



Figure 2. 68%, 90% and 99% limits (blue, red and green line 
respectively) on the unresolved line normalization in the 0.5-5 
keV energy interval. 



ral filtering, event screening, etc.), can be found in Owen & 
|Warwick~||2010| We wish to emphasize here that the data 
were spatially filtered in order to remove the contamination 

erg 



due to bright sources down to a luminosity of 2 x 10 35 



s~ . A residual contamination is estimated by the authors 
which, as done in the referred paper, is taken into account 
in the following spectral analysis. 

Appropriate Auxiliary Response File and Response Ma- 
trix File for the source-filtered region were created, while 
background files were produced from "blank-sky" fields ex- 
tracted from a region of sky close to M33, and "filter-wheel 
closed" data to produce a background spectrum (again see 
details in Owen fc Warwick |2010[ ) . The final spectra for the 
three observations are shown in Figure [l] 

We first verified that we were able to reproduce the 



Owen & Warwick 2010 results, fitting (using the Cash 



statistics) the data with XSPEC 12.6 softwartj^] using a 
multi-component model consisting in two thermal plasma 
subject to Galactic absorption of Nh = 7.5x 10 20 cm -2 , con- 
strained to solar abundance, with temperatures of kT = 0.2 
keV and kT — 0.6 keV. We also tested the use of models 
with subsolar (Zq — 0.1) metallicity but the results are con- 
sistent with those assuming solar metallicity. In addition a 
power law component was added to account for the 9% flux 
from the bright sources leaking outside the masked regions; 
the bright ULX in the galaxy center, which can be fitted 
with a more complex spectrum (power law, disk black body, 
and intrinsic absorption) was treated separately, assuming 
that 3% of its flux spills into our source-filtered region, as 
done in Owen & Warwick 2010 We then added an unre- 
solved line to our spectral model, and fitted iteratively the 
whole spectrum in the range 0.5-5 keV allowing the line flux, 
and the thermal components parameters (temperature and 
normalization), free to vary. 

We computed the 68%, 90% and 99% confidence levels 
for two interesting parameters (line flux and central energy), 
shown in Figure [2j sampling our energy range in steps of 20 
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Figure 3. Ratio between the flux upper limits derived with the 
old (SAS v8) and new (SAS vll) calibration. Calibration issues 
introduce at most a factor of two uncertainty in the result. 



eV, in order to fully exploit the EPIC spectral resolution. 
We note that several residual features are visible in the fitted 
spectrum, such as the emission lines around 1.5 keV due to 
A1-Kq line. Furthermore, the spectral model we use does 
not account for unresolved point sources below the 2 x 10 35 
erg s _1 threshold; this contribution is estimated by | Owen fc| 



Warwick] ( |2010[ ) to be < 1% of the bright source emission. In 
order to minimize the uncertainties involved in refining the 
background spectrum, or in extrapolating the point source 
luminosity function to lower fluxes, we prefer to leave these 
contribution in the data, thus obtaining conservative lower 
limits on the DM line flux. 

For consistency, as well as because of the very detailed 



checks performed in the Owen & Warwick (20101 work, we 



prefer to use here the original data. Nonetheless, we checked 
the robustness of our results with respect to an update of 
the calibration software repeating the analysis -initially per- 
formed with SAS v8- with the new SAS vll. The ratio be- 
tween the flux upper limits derived with the old and new 
calibration is plotted in Figure [3] and shows that calibra- 
tion issues introduce an uncertainty, at most, of a factor of 
two in the result. We reasonably expect that further releases 
of the calibration software will not change the result more 
than this. 



3 X-RAY FLUX FROM STERILE NEUTRINO 
DARK MATTER 

We use the upper limits on the X-ray line flux derived in 
the previous section to constrain the possible contribution 
to M33 X-ray emission due to a radiative decay of sterile 
neutrino DM. This is done assuming that this additional 



contribution to the flux cannot exceed the upper limit de- 
duced in the previuos section. 

Majorana sterile neutrinos may decay to a photon and 
an active neutrino ( Pal fc Wolfenstein 19821. The decay 
width of the process i^] 



F(u s JV a ) = 1.36 X 10 



-32 sin 2 (2(9) / m B \ ' 

10- 10 Uevy 



where 6 is the neutrino mixing angle and m„ -the double 
of the emitted photon energy- is the mass of the sterile 
neutrino. 

The resulting flux at the Solar System position is $ = 
F SI S/47rm s , in which f2 and S represent the angular size 
and the column density of the emitting region, respectively. 
In order to evaluate the DM column density, defined as the 
integral of the mass density along the line of sight, we need 
to parametrize the matter content of M33. Since there is no 
universal consensus on the DM halo profile in M33, we tried 
both the cored and spiked profiles that are consistent with 
the galaxy rotation curve ( jCorbelli 12003) . T he cored case is 
well represented by a Burkert density profile ( Burkert |1995 \ 

p(r)=p (l + x)- 1 (l + x 2 y\ 

while the spiked profile is well described by a NFW density 



distribution (Navarro et al. 
p(r) = p x' 



1997) 



-\l + x) 



with x = r/ro. We consider two extremal cases among the 
various fits obtained in Corbclli 2003, taking the 3a lower 
limit for the DM density in the Burkert case, and the 3cr 
upper one in the NFW case. This way we are taking into 
account corresponding 3<r variation of the M/L ratio in the 
inner part of M33 ( jCorbelli |[2003| ). Their defining parame- 
ters are: 



model 


r (kpc) 


po (GeVc- 2 cm- 3 ) 


NFW 


35 


0.0765 


Burkert 


12 


0.336 



Profiles with a central slope steeper than r 1 like the Moore 
profil e ( jMoore et aT~||1999 l are instead excluded (ICorbelli 
2003k. 



Figure [4] shows the DM column density for the various 
parametrization of M33 considered here: filled triangles in 
the NFW case, empty triangles for a Burkert profile. We 
assumed an infinite halo in calculating the mean surface 
density, because we checked that, repeting the calculation 
for a 17 kpc radius halo (the minimum possible radius, cor- 
responding to the furthest observation in |Corbelli~ 2003 \ 
the result would change by a factor ~ 1.3. In each case we 
considered separately the inner region (in) of the galaxy - 
corresponding to the central 3.5 kpc bulge- and the outer 
part (out) of it -defined as the cylindrical volume corre- 
sponding to the same angular region and from which the 
central bulge has been subtracted. In addition we also con- 
sidered the sum of both contributions (tot). In all the cases 
we took into account only the flux coming from the ellipti- 



cal region analyzed in Owen & Warwick (2010 1. The values 



2 For Dirac sterile neutrinos the decay width of is half the one 
assumed here I Pal fc Wolfenstein [1982 I. 
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Figure 4. Column density for both mass paramctrization of M33 
considered here compared to the ones assumed in previous similar 
analysis. The angle gives the angular distance from the galactic 
center (126.1° in the case of M33). The solid line represents the 
column density of the Milky Way (see text for more details). 



we obtain for the column density range from 40.4 to 266 
Mq/pc 2 . This precaution addresses the concerns raised by 
|Kusenko fc Loewenstein"] ( |2010[ ) about the uncertainties af- 
fecting the DM profile in the central region of the galaxy 
(see next section). A comparison is made with the column 
density of the Milky Way and the other systems considered 
in |Loewenstein fc Kusenko |201(j] and|Boyarsky e t al. |2010| 
In evaluating the total flux at the Solar System position we 
added up the contribution due to the Milky Way along the 
M33 line of sight, evaluated according to the DM density 
parametrization given in Strigari e t al. |[2~008| 

Figure [5] shows the regions of the parameter space 
(sin 2 20, m s ) excluded by the X-ray observations of M33. 
The main assumption here is that the DM halo is entirely 
made up by sterile neutrinos. Limits are shown at 1, 2, and 3 
a confidence level. We show the most and the least stringent 
cases, among the six we considered here: NFW tot (dotted 
lines) and Bur. in (solid lines), respectively. The yellow line 
gives the models in which the cosmological amount of DM 
is entirely made up by sterile neutrinos. The yellow region is 
therefore excluded by over-closure. The star represents the 



model of Loewenstein & Kusenko (20101 which is found to 



be inconsistent with the observed X-ray emission from M33. 



4 DISCUSSION AND CONCLUSIONS 

Let us now summarise the assumptions under which our 
results were obtained and the uncertainty affecting them. As 
pointed out in Kusenko & Loewenstein 2010, the interaction 
of DM with baryons could have the effect of expelling the 
DM out of the central part of a spiral galaxy ( Kl ypin "et| 
al. |200l" I . DM could therefore provide only a sub-dominant 
contribution to the total amount of mass in the central part 
of M33. This is the reason why we separated the emission 
along the line of sight in two contributions, the one produced 
in the innermost region of M33, and the contribution coming 
from the two external cylinders. Interestingly (see Figure Hk, 
these external (out) contributions -whose mass distribution 



is determined on a more stable basis- are the ones giving 
the most stringent bounds. This makes our results stable 
against the main criticism raised by Kusenko & Loewenstein] 
[] ( |2010[ ) . We stress that our results change by less than one 
order of magnitude because of the uncertainty affecting the 
DM density distribution, that could be either cored or spiked 
at the center of M33. 

An additional uncertainty factor is represented by the 
effect of photometric absorption. This will reduce the de- 
tected X-ray flux, thus resulting in an underestimate of the 
M33 DM emission. We note however that: i) galactic ab- 
sorption has already been taken into account in the spectral 
fitting described in Sec. 2; ii) photometric absorption is rele- 
vant only in the soft X-ray band below 2 keV, and thus the 
limits in the hard energy range are essentially unaffected; 
iii) additional intrinsic M33 absorption may be present but 
likely with a low covering factor. The population of M33 X- 



ray sources for instance (see e.g. discussion in Foschini et al 



2004 and Grimm ct al. 200T] does not reveal widespread 
absorption, and anyway with nu less than a few 10 21 cm -2 
which results in a factor of ~ 2 reduced flux at 1 keV. Nev- 
ertheless, to evaluate the effect that absorption might have 
on our estimates, we could simply divide by ~ 2 the emis- 
sion coming from the back cilinder. This would result in a 
25% reduction of the total flux (out cases), which does not 
affect our main results. We also need to point out that in our 
analysis we did not try to fine tune the background removal, 
nor we did model the contribution of the unresolved X-ray 
binary population, in order to avoid biases due to the poor 
knowledge of these contributions to the overall X-ray emis- 
sion. This choice results in conservative upper limits on the 
M33 DM content, since the estimated emission is actually 
overestimated. 

In summary, we make use of archival XMM-Newton ob- 
servation of the diffuse and unresolved components emission 
of the inner disc of M33 ( |Owen fc Warwick |[2010] | , finding 
no evidence for line emission compatible with sterile neu- 
trino radiative decay. We accordingly set bounds on sterile 
neutrino parameter space under the assumption that they 
account for the entire DM content of M33. In particular we 
have found no evidence for the emission by 5 keV neutrinos 



deduced by Loewenstein & Kusenko (2010) 
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